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Abstract — Despite the occurrence of elegant algorithms for 
solving complex problem, exhaustive search has retained its 
significance since many real-life problems exhibit no regular 
structure and exhaustive search is the only possible solution. 
The advent of high-performance computing either via multicore 
processors or distributed processors emphasizes the possibility 
for exhaustive search by multiple search agents. Here we analyse 
the performance of exhaustive search when it is conducted 
by multiple search agents. Several strategies for cooperation 
between the search agents are evaluated. We discover that the 
performance of the search improves with the increase in the level 
of cooperation. Same search performance can be achieved with 
homogeneous and heterogeneous search agents provided that the 
length of subregions allocated to individual search regions follow 
the differences in the speeds of heterogeneous search agents. 

Index Terms — Parallel algorithms, Search methods. 



I. Introduction 

Exhaustive search consists of systematically enumerating 
all possible candidates for the solution and checking 
whether each candidate satisfies the problem's statement. The 
number of candidate solutions to consider grows very rapidly 
with problem size, causing lengthy or even infeasible searches. 

The continuing increase in computing power and memory 
sizes, and the advent of multicore processors and parallel and 
distributed programming (TJ, 12 (combined effort of huge 
number of computers connected on the internet) increases 
the feasibility of exhaustive search and has revived interest 
in brute-force techniques for a good reason. Many real-life 
problems reveal no regular structures to be exploited in the 
search for solutions, and this leaves exhaustive search as the 
only possible approach. We can always try to design elegant 
and optimal algorithms in a quest for order of magnitude of 
performance improvements. However, mathematical creativity 
is not guaranteed to give success for real-life problems. As 
pointed out too many times in the past, it is more likely 
that order of magnitude improvements can be achieved due 
to program optimization, and to the clever use of limited 
computational resources in an exhaustive search. 

There are many problems where the exhaustive search 
is the best solution. Radio SETI (Search for Extraterres- 
trial Intelligence) uses radio telescopes to listen for narrow- 
bandwidth radio signals from space, which are not known 
to occur naturally, so a detection would provide evidence 
of extraterrestrial technology. Radio SETI has an insatiable 
appetite for computing power to cover greater frequency 
ranges with more sensitivity. Early radio SETI projects have 
used special-purpose supercomputers, located at the telescope, 



to do the bulk of the data analysis. SETI@home is doing 
radio SETI using a virtual supercomputer composed of large 
numbers of Internet-connected computers O . It was originally 
launched in May 1999. As of 2008, five million people 
in 226 countries have participated volunteered their PCs to 
analyse data. Combined, their PCs form EarthaAZs second 
most powerful supercomputer, averaging 482 TeraFLOPs and 
contributing over two million years of CPU time |4). 

DES (Data Encryption Standard) Q is a secret-key block 
cipher that was selected by the National Bureau of Standards 
as an official Federal Information Processing Standard (FIPS) 
for the United States in 1976. Its widespread international 
use gave rise to extensive cryptanalysis of DES. None of the 
proposed cryptanalytic attacks was practically feasible, and 
eventually it was broken using exhaustive search of its 56-bit 
key region. In 1998, Electronic Frontier Foundation [6] built 
a custom-made machine named "Deep Crack" for the cost of 
US$250,000, and managed to decrypt the DES Challenge II-2 
test message after 56 hours of work. COPACOBANA machine 
built in 2006 Q is the only other confirmed DES cracker, but 
unlike Deep Crack, COPACOBANA consist of commercially 
available, reconfigurable integrated circuits. COPACOBANA's 
successor named RIVYERA reduced the time to recover a 
DES key to the average of a day (8). 

In this paper we evaluate the performance of several 
strategies for exhaustive search depending on the following 
parameters: 

• Homogeneous or heterogeneous search agents according 
to their speed. What is the impact of the differences 
in speeds of search agents on the overall speed of the 
search? 

• Length of allocated search regions. What is the impact 
of the division of the search region on the speed of 
cooperative search? 

We will also demonstrate the optimum division of the search 
region. 

Here is the overview of the paper. In Section [TT] we present 
theoretical results on the average search time for exhaustive 
search. Section [III] we analyse exhaustive search with homo- 
geneous cooperating agents, while in Section [IV] we analyse 
cooperation of heterogeneous agents. Section fVl concludes the 
paper and gives directions for future research. 

II. Exhaustive search methods 

First we give a formal definition of exhaustive search. 
Consider function F : X — > [0, 1] where X is the discrete 



domain of F. X is also called the search region. Assume 
that there is only one point x s £ X such that F(x s ) = 1. 
Otherwise, F(x) = if x £ X and x ^ x s . We call 
point x s solution of function F. If function F exhibits no 
regular structure, then finding the solution x s by exhaustively 
searching the domain R can be the only option. Exhaustive 
search attempts to find the point x s by repeatedly calculating 
F(x) for all points X £ R until x s is found such that 
F(x s ) = 1. We assume that the region is cooperatively 
searched by m agents ai, a,2, • • • , a m . We are interested in 
the impact of the cooperation between multiple search agents 
on the performance of the search. Search performance is 
measured bu the average time required to find the only solution 

J g . 

For sake of simplicity and clarity, we analyse exhaustive 
search of one-dimensional region X = [1,L]. However, results 
are equaly valid for multi-dimensional regions since they can 
be easily converted into a one-dimensional region. Each agent 
a.i is allocated a subregion Xj C X which it searches with 
speed Vi. An obvious requirement is that Ui=i™ -^t = X. Each 
agent starts its search from the starting point of the allocated 
subregion. If an agent reaches the end of the search region, 
then it continues from the beginning of the search region. If 
m-th agent reaches point L, then it continues the search from 
point 1. This is illustrated in Fig. Q] 
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Fig. 1. Exhaustive search with multiple search agents. 



The process of allocation of search subregions to search 
agents depends on the type of communication between search 
agents and the central server, whether the central server knows 
the number m of search agents and the speed of each agent, 
and whether the number of search agents is determined before 
the search starts or additional agents can join the search at a 
later time. One possibility is that the central server can allocate 
a subregion to each search agent depending on the number 
m of search agents and the speed of each agent. Another 
possibility is that each agent chooses its own starting point 
for the search irrespective of the starting points of the other 
search agents. Search continues until a solution is found, or a 
"stop" command is received from the central server. 

We assume uniform probability distribution of the solution 
in the search region X, that is, each point in the search region 
X is equally probable to be the solution x s . Average search 
time is calculated using the following formula 
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(1) 



where (i) Pr{x s £ i?/} is the probability that the solution is 
in region R[ with length I; (ii) p v (v) is the probability that the 
speed of a search agent is v; and (iii) l/(2v) is the average 



time to find the solution in a region of length I by agent with 
speed v. 

In Eq. (Q]i it is assumed that the two random variables are 
mutually independent. If not, then in Eq. (Q]i one needs to use 
joint probability distribution function 
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If the search region with length L is divided into m 
subregions, and length I of subregions is a random variable 
with probability distribution function pi(l), then one can easily 
show that 

mpi(l)l 



Pr{x s £ R t } 



(3) 



where mpi(l) is the number of regions with length I, mpi(l)l 
is the total number of points belonging to regions with length 
I, and the ratio mpi(l)l/L is the probability that the solution 
x s is in a region with length I. Substituting Eq. (0) in Eq. ©, 
one obtains 
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If v and I are mutually independent random variables, then 
Eq. transform into 



III. Exhaustive search with homogeneous agents 

In this section we additionally assume that all search agents 
are homogeneous, that is, all search agents are with same 
f 1 for v = V 

speed V andp„(w) = < . Here we compare three 

I otherwise 

methods for division or the search region, which are explained 
next. 

A. Equal subregions 

This can occur if the number of search agents m is known in 
advance, and two directional communications exists between 
central server and search agents. Thus a search subregion is 
allocated and communicated to each search agent in advance. 
In this case / 

1 for I = i 

m 

otherwise 



pi(i) = 



(6) 



Thus Eq. <(5j yields 
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B. Semi-equal subregions 

Search subregion is allocated and communicated to each 
search agent by the central server as it registers. Then search 
subregion allocated to a newly registered agent depends on the 
current number of search agents. This case can occur if the 
number of search agents is not known in advance, and two 
directional communications exists between central server and 
search agents. New search agents can register at run time. Let 



[1, L] denote the entire search region. Then the first agent will 
start its search from point 1 . If there is only one search agent 
in the system, then it searches the whole range. If there is a 
second search agent in the system, then its searching subregion 
is the second half of the whole region, that is, the second 
agent will start the search from point l + L/2, If the searched 
solution is located in the second half of the range, then the 
second search agent will find it faster than the first search 
agent. Third agent will start from point l + L/4, and the next 
agent will start from point 1 + 3L/4. The next agents will 
start the search from points 1 + L/8, 1 + 3L/8, 1 + 5L/8, 
1 + 7L/8 etc. Thus, the searching subregions are shrinking as 
the number of search agents grows. Each time a new search 
agent joins the search, it is given to search the second half 
of the currently largest search subregion. In the semi-equal 
subregions method, for m = 2°, 2 1 , 2 2 , 2 3 , . . . , p;(Z) and E{t) 
are again given by Eq. © and Eq. Q. Hence the semi-equal 
subregions method is identical to the equal subregions method. 
Otherwise, for 2" < m < 2" +1 , it is easy to show that the 
probability that a search subregion is with length I is given by 



2 „H 



w(0 







for I = £ 
for I = ^ 
otherwise 



(8) 



The average search time calculated using Eq. ([8]l is higher 
than the average search time for the equal subregions method 
obtained from Eq. ||7). 

C. Random subregions 

Each search agent starts from a randomly chosen starting 
point, thus randomly choosing its search subregion. Therefore, 
the size of the subregion searched by an agent can vary 
between and the size L of the entire search region. This 
kind of search can occur if, for example, one directional 
communication exists from search agents to central server. 
Number of search agents is not known in advance, and 
new search agents can join at run time. The one directional 
communication is used by a search agent to communicate to 
the central server when the solution is found. 

For given L, m and when agents randomly choose start- 
ing points independently from each other, it is possible to 
analytically calculate pi(l). However, our main interest in 
this paper is the average search time and comparison with 
the other two methods Equal subregions and Semi-equal 
subregions. Therefore, we have taken the simpler approach 
and have numerically calculated the probability function pi(l) 
for L = 1000 and m = 2, 3, 4, . . . , 32 using the Monte-Carlo 
method [9 1. Figure|2]depicts the calculated probability function 
Pl(l). For m = 2, pi(l) is uniformly distributed probability 
function in the region [0,999]. For m = 3, pi(l) is linearly 
decreasing function in the region [0,999]. For m > 3, pt(V) 
is monotonically decreasing polynomial function in the region 
[0,999]. 

D. Comparison 

Figure |3] gives the average search time (y-axis) for the three 
methods depending on the number m of homogeneous search 
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Fig. 2. Numerically calculated probability function p; (I) for L = 1000 and 
m = 2,5, 10,20,30. 
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Fig. 3. Comparison of performances of the three searching methods. 

As expected Equal subregions method produces the best 
performance, i.e., the shortest average search time. Semi- 
equal subregions method produces performances which are 
close to the ones produced by the Equal subregions method. 
Random subregions method results in significantly higher 
average search times. For example, 10 search agents using 
the equal subregions method will produce same performance 
as 19 search agents using the random subregions method, 
and 16 search agents using the equal subregions method will 
produce same performance as 31 search agents using the 
random subregions method. 

IV. Exhaustive search with heterogeneous agents 

In this section we consider the performance of exhaustive 
search when the cooperating agents are heterogeneous, that 
is, agents search the search region X with different speed. 
Average search time is calculated using Eq. (01 and Eq. ©. 
An interesting consequence from Eq. (0 is following: if v and 
I are mutually independent random variables, then average 
search time depends on E(l/v). In other words, heteroge- 
neous agents will provide the same search performance as 
homogeneous agents whose search speed is V — E(l/v). 



In our experiments, searching agents have speed which is 
a random variable with the following probability function: 
P„(0.5) = 0.3, P„(1.0) = 0.3 and P„(1.375) = 0.4. One 
can easily calculate E(l/v) = 1.19091. 

We analyse the performance of exhaustive search with co- 
operating heterogeneous agents for the following three search 
strategies. 

A. Strategy 1: One-directional search 

Each agent is randomly allocated a starting point for the 
search. Then the agent searches the region from its starting 
point untill the starting point for the next agent. Starting point 
is a uniformly distributed random variable in the region [1,L]. 
Length of allocated search subregions is a random variable 
whose probability function is depicted in Fig. |2] Then each 
agent is searching its own subregion (see Fig. [TJ. It is possible 
that agent at finishes the search of its region before the 
solution x s is found by any of the agents. Then agent a, 
continues with the search of the subregion until a "stop" 
command is received by the central server. Let v m i n and v max 
denote the minimum and maximum speed of search agents, 
and l m i n and l max denote the minimum and maximum length 
of search subregions. Then provided that 

Vmax ^ lmin ~t~ Imax 
^min lmax 

agent a.j+1 will finish searching subregion Xj+i before agent 
di will finish searching both regions and JQ+i for all i = 
1,2, ... ,m. In other words, each agent is responsible for the 
search of its allocated subregion and will receive no help from 
other agents in the search of its subregion. 

B. Strategy 2: Two-directional search 

Similar to the previous case, each agent is randomly al- 
located a search subregion by means of an assigned starting 
point. However, neighbouring agents can help each other in 
the following manner: each agent a% conducts the search of 
its subregion Xi in two directions: to the left and to the right 
from the assigned starting point (see Fig. @}. If the speed of 
search for agent a* is Vi, then the search to the left side 
is conducted with speed v±/2. Same speed applies for the 
search to the right side. Searching in both directions improves 
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Fig. 4. Exhaustive search with cooperative agents in two directions. 

the cooperation between neighbouring agents. Assume that the 
distance between the starting points for agents cii and a 2 is 
di 2, and the distance between the starting points for agents 
a 2 and a 3 is <i 2 3 . If a 2 is faster than a\, then a 2 will search 
larger portion of d\ t2 than o-i- Similarly, if 0,2 is slower than 



a 3 , then a 2 will search smaller portion of e? 23 than a 3 . In 
other words, agents help their slower neighbours, and get help 
from faster neighbours. It is rather straightforward to show 
that, for example, agent a 2 searches a subregion with length 
h = hi + where l 21 = ±^v 2 and l 23 = ±^v 2 . Time 
spent by agent a 2 to search the subregion X 2 is given by the 
following equation: t 2 = (7 2 i + hz)/v2. 

Figure [5] gives the average search time (y-axis) for the three 
methods depending on the number of heterogeneous search 
agents (x-axis). It is obvious that the search strategy with 
two neighbours helping each other (two-directional search) 
significantly reduces the search time compared to the one- 
directional search strategy where the neighbouring agents do 
not help each other. An obvious question arises: can the search 
time be further reduced if agents help each other in groups 
of three, that is, agents in groups of three jointly search 
an allocated subregion? Defining the strategy for cooperation 
amongst groups of three neighbours and for joint search of an 
allocated subregion is beyond the scope of this paper. Here 
we are interested only in the search performance. Figure [5] 
confirms our intuitive expectations: joint search of allocated 
regions in groups of three agents further reduces the average 
search time. Further increasing the number of agents that 
jointly search a subregion to four additionally reduces the 
search time, as shown in Fig. [5] too. 
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Fig. 5. Average search time for cooperating heterogeneous search agents, a) 
random; b) joint search in groups of two agents; c) joint search in groups of 
three agents; d) joint search in groups of four agents; e) length of subregion 
is determined by agent's speed. 

Figure [6] gives a possible strategy for joint search by groups 
of four search agents in a two-dimensional search space X. 
Each agent should search in four directions simultaneously, 
that is, in two directions for each of the two dimensions. If 
the speed of agent en is Vi, then the search speed in each 
of the four directions will be i>j/4. If a subregion ^12.3 4 
is jointly searched by four agents 0,1,0,2,0,3,0,4 with search 
speeds V\, v 2 , 1)3, V4, respectively, then agent a\ will search 
v i/{ v i + v 2 + v 3 + V4) parts from region -X"i,2,3,4> agent a 2 
will search v 2 /(vi+v 2 -\-V3+V4) parts from region ^1,2,3 4 etc - 
Subregion X 22 searched by agent a 22 is shown as a shaded 
square in Fig. [6] Careful examination of the boundaries of 



the subregion X 22 reveals that v 2 2 < w 23 ,w 2 2 > "32,^22 < 

Vl2,V 2 2 = V 2 1- 
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Fig. 6. Exhaustive search in a two-dimensional search region and cooperation 
with four neighbours. 

Figure [6] also gives a hint on a possible strategy for coop- 
eration between 2 N neighbouring agents, for TV = 2, 3, If 

the search region is TV-dimensional, then each agent can be 
assigned a random starting point and the agent can start the 
search in 2N directions - two directions for each dimension. 
If the dimension of the search region X is less than N, e.g. 1, 
then one can transform X into a TV-dimensional search region. 

C. Strategy 3: Subregion 's length is proportional to agent's 
speed 

Each search agent is allocated search subregion whose 
length is proportional to agent's speed. Faster agents get larger 
subregions; slower agents get smaller subregions. The size k 
of the subregion Xi allocated to agent a t is determined by the 
following equation: 
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(10) 



Vi + V 2 + ■ ■ ■ + V m 

Each agent will finish the search of its allocated subregion at 
the same time 

L 



k 

U = - 

Vi 



V1+V2 + 

and the average search time is given by 
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E(t) 



(11) 



(12) 
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This strategy is the optimum strategy for division of the search 
region for heterogeneous agents. 

D. Comparison 

As we see from Fig. [5] in the case of heterogeneous agents, 
average search time decreases and the search performance 
improves as the cooperation between agents grows. One- 
directional search is the worst search strategy: cooperation is 
reduced to the division of the search region between the search 
agents. If faster agents are enabled to help slower agents e.g. 
by searching the allocated one-dimensional subregion in two 
directions, then the average search time reduces dramatically. 
If the number of agents n that jointly search a subregion 
grows, then the averages search time further reduces. For 



example, one-directional search with m = 23 agents, two- 
directional search with m = 14 agents for groups of two 
agents, two-directional search with m = 12 agents for groups 
of three agents, two-directional search with in = 11 agents 
for groups of four agents, and optimum search with m = 10 
agents produce simillar average search times. For n — m, two- 
directional search strategy converges to the optimum strategy 
and each agent searches a subregion whose size is proportional 
to the agent's speed. 

We also note that the optimum strategy (Subregion 's length 
is proportional to agent's speed) for heterogeneous agents 
produces same search performance as the Equal subregions 
method for homogeneous agents, if E(v) for the heteroge- 
neous agents is equal to the search speed V of homogeneous 
agents. 

V. Conclusion 

We have analysed the performance of exhaustive search by 
cooperative search agents. Both homogeneous and heteroge- 
neous agents are analysed. Performance of exhaustive search 
by cooperative search agents improves and average search 
time decreases as the level of cooperation increases. Optimum 
performance is achieved if the central server knows the number 
and speed of search agents, and then each agent is allocated 
a search subregion with length proportional to agent's speed. 
If the search region is with high-dimension TV close to the 
number of agents m and two-directional search strategy is 
employed, then search performances close to the optimum 
can be achieved. Results given in this paper are presented for 
exhaustive search, but they are equally valid for other search 
methods where multiple search agents cooperate. 
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